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Abstract 

Diverging junctions are important network bottlenecks, and a better understanding of di- 
verging traffic dynamics has both theoretical and practical implications. In this paper, we first 
introduce a continuous multi-commodity kinematic wave model of diverging traffic and then 
present a new framework for constructing kinematic wave solutions to its Riemann problem 
with jump initial conditions. In supply-demand space, the solutions on a link consist of an inte- 
rior state and a stationary state, subject to admissible conditions such that there are no positive 
and negative kinematic waves on the upstream and downstream links respectively. In addi- 
tion, the solutions have to satisfy entropy conditions consistent with various discrete diverge 
models. In the proposed analytical framework, kinematic waves on each link can be uniquely 
determined by the stationary and initial conditions, and we prove that the stationary states and 
boundary fluxes exist and are unique for the Riemann problem of diverge models when all or 
partial of vehicles have predefined routes. We show that the two diverge models by Lebacque 
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and Daganzo are asymptotically equivalent. We also prove that the supply-proportional and 
priority-based diverge models are locally optimal evacuation strategies. With numerical exam- 
ples, we demonstrate the validity of the analytical solutions of interior states, stationary states, 
and corresponding kinematic waves. This study presents a unified framework for analyzing 
traffic dynamics arising in diverging traffic and could be helpful for developing emergency 
evacuation strategies. 

Key words: Kinematic wave models, diverging traffic, Riemann problem, supply-demand space, 
stationary states, interior states, boundary fluxes, turning proportions, First-In-First-Out, evacua- 
tion strategies 



1 Introduction 

Essential to effective and efficient transportation control, management, and planning is a better 
understanding of the evolution of traffic dynamics on a road netwo rk, i.e., the formation, propa 



gation, and di ssipation of traffic congestion. The seminal work by (Light hill and WhithamL 



Richards, 



1955 : 



1956) (LWR) describes traffic dynamics with kinematic waves, including shock and rar- 
efaction waves, in density (p), speed (v), and flux (q). Based on a continuous version of traffic 
conservation, ^£ + |^ = o, and an assumption of a speed-density relationship, v = V (p), the LWR 
model can be written as 

|p + |-pV<p)=0, (1) 

which is for a homogeneous road link with time and location independent traffic characteristics, 
such as free flow speed, jam density, capacity, and so on. In general, V(p) is a non-increasing 
function, and vy = V(0) is the free flow speed. In addition, q = Q(p) = pV(p) is unimodal with 
capacity C = Q(p c ), where p c is the critical density. Traffic states with density higher than p c are 
congested or over-critical, and those with lower density are free flowing or under-critical. Here we 
denote the jam density by pj, and p G [0,py]. 



In a road network, however, more important and interesting are the formation, propagation, and 
dissipation of traffi c queues caused by n etwork bottlenecks, including merges, diverges, and other 



network junctions (|Daganzo et al 



1999). But compared with numerous studies on the LWR mode l 



and higher-order models of traffic flow on a road link (Fe deral Highway Administr atio 



ion! 



2004), 



studie s on traffic dynamics at merging, diverging, and other junctions are scarce. In (|Fazio et al. 



990 ), behavioral mo dels were proposed to capture individual vehicles' diverging maneuvers. In 



(|Papageor giou. 



1990), diverging flows of vehicles on a path are determined by pre-defined splitting 

"J 



rates. In (|Daganzol ll995). the First-In-First-Out (FIFO) principle was explicitly introduced so that 
diverging flows are proportional to turning proportions, which can be time-dependent. But it was 
noted th at the FIFO principle could be vio lated when one downstream branch is heavily congested. 



In (Liu et al 



19961 : iNgoduv et al. 



2006), diverging traffic was considered in a so-called friction 



term of a higher-order model, where diverging flow to an off-ramp is determined by e xpected 
diverging flow and the congestion level of the off -ramp. In (jMunoz and DaganzoL |2002|) . it was 



shown that First-In-First-Out (FIFO) blockage caused by one congested downstream branch could 
significantly reduce the discharging flow -rate of the whole diverge, and vehicles may not f ollow the 



FIFO 



1997 



principle strict 



Daganzo et al. . 



. Diverging traffic w ith t wo or more veh icles have been studied in (Daganzo 



19971 : 



Newell . 



1999). In(Cassidv, 



2003), metering strategies were discussed 



1999J), different network bottlenecks can 



for diverging junctions. As pointed out in (|Daganzol . 
induce different traffic behavior; at diverging junctions, which are different from merging and 
other junctions, not only capacities of all branches but also the combinations of diverging vehicles 
on the upstream branch could determine the formation and dissipation of queues. In addition, a 
better understandin g of diverging traffic flow could also lead to more efficient evacuation strategies 



(Sheffietal 



19820 . In this study, we are interested in traffic dynamics arising from diverging 



junctions for one type of vehicles within the framework of the LWR model. 

Considering the analytical power and simplicity of the LWR model, many researchers have 
attempted to study traffic dynamics arising in general transportation networks in the framework of 
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kinematic wave models. In one line, 



Daganzol(| 19951) and 



Lebacqud (|1996l) extended the Godunov 



discrete form of the LWR model for computing traffic flows through diverging, diverging, and 
general junctions. Hereafter we call such models as Cell Transmission Models (CTM). In CTM, 
so-called traffic demand and supply functions are introduced, and boundary fluxes through vari- 
ous types of junctions can be written as functions of upstream demands and downstream supplies. 
In CTM, various physically meaningful rules can be used to compute boun dary fluxes, such as 



the First-In-Fi rst-Out diverging princ iple ( 



ing principle (|Jin and Zhang . 



PapageorgiouL 



1990; 



Daganzo . 



1995j) and the fair merg- 



2003bl) . CTM are discrete in nature and only suitable for numerical 



simulations. Thus they do not provide any analyt ical insights on traffic dyna mics at a network 



intersection as the LWR model. In another line, 



Hol den and Risebrol (| 19951) and 



Coclite et al. 



(I2005|) attempted to solve a Riemann problem of an intersection with m upstream links and n 
downstream links. In both of the analytical studies, all links are homogeneous and have the same 
speed-density relations, and traffic dynamics on each link are described by the LWR model. In 



(|Holden and Risebrol . 



1995), the Riemann problem with jump initial conditions is solved by in- 
troducing an entrop y condition that maximizes an objective function of all boundary fluxes. In 
(|Coclite et all 120051) . the Riemann problem is solved to maximize total flux with turning propor- 
tions. Both studies were able to describe basic waves arising from a network intersection but 
also subject to significant shortcomi ngs: (i) All links are as sumed to have the same fundamen- 



tal diagram in both studies; (ii) In (jHolden and Risebro . 



19951) . vehicles can travel to an arbitrary 



downstream link, and the entropy co nditions used are pragmatic and lack of physical interpreta- 



tions; and (iii) In (ICoclite et all 120051) . results are only valid for restricted turning proportions and 



junctions with no fewer downstream links; i.e., n > m. In addition, neither of these studies present 
a unified continuous model of network vehicular traffic. 



As in (Holden and Risebro 



1995 



Coclite et al. 



2005J), in this study we attempt to analytically 



obtain kinematic wave solutions of traffic dynamics arising at a diverging junction. However, our 
study does not bear the same limitations as in these studies: all links can be mainline freeways or 
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off-ramps with the same or different characteristics, and our solutions are physically meaningful 
an d consistent with the discrete supp ly-demand models of diverging traffic, e.g. those proposed 



in (|Daganzo 



1995 



Lebacqud . 



19961) . We first present a continuous kinematic wave model of 
multi-commodity diverging traffic flow based on the conservation of commodity traffic. Following 
the new framework used to solve Riemann problems for in homo geneous LWR model at a linear 



junction (|Jin et al. 



2010) . we present a new framework 



2009) and for merging traffic flow (Jin, 
for solving the Riemann problem for diverge models. In the Riemann solutions, there can be a 
stationary state and an interior state for each branch. Here stationary states are the self-similar 
states at the boundary. That is, in the Riemann solutions, stationary states prevail all links after 
a long time. In contrast, interior states do not take any space in the continuous solution and only 



1984). We introduce a 



show up in one cell in the numerical solutions as observed in (|van Leei . 
so-called supply-dem and diagram and discuss the problem in supply-demand space, rather than 



in p — q space as in (IHolden and Risebrd . 1 19951 : ICoclite et al. 



20051) . After deriving admissible 



solutions for upstream and downstream stationary and interior states, we introduce an entropy 
condition based on various diverge models. We then prove that stationary states and boundary 
fluxes are unique for given upstream demand and downstream supplies (but interior states may 
not). Then, kinematic waves on a link are determined by the corresponding LWR model with the 
stationary state and the initial state. In a sense, kinematic waves of the Riemann problem can be 



rules in ( 
Diffe 


Daganzo. 


1995; 


Lebacaue, 


199( 


). 


rent from 


(Holden and Risebro. 


1995; 


Coclite et al. 



2005), where the Riemann solutions 



only comprise of initial and stationary states, here we have additional interior states. Interior states 



were obse rved when the i nhomogeneous 



ring road (IJin and ZhangL 



2003a; 



Jin et al. 



WR m odel was used to simulate traffic dynamics on a 



2009). Although interior states are not directly related 



to kinematic waves on all links, they are used in the entropy condition and therefore essential to 
picking out unique physical solutions. As we can see later, interior states are essential to construct 
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kinematic wave solutions for different diverge models. 

The rest of the paper is organized as follows. In Section 2, we introduce a continuous multi- 
commodity kinematic wave model of diverging traffic. In Section 3, we introduce a new framework 
for solving the kinematic waves of the Riemann problem with jump initial conditions in supply- 
demand space. In particular, we derive traffic conservation conditions, admissible conditions of 
stationary and interior states, and additional entropy conditions based on various discrete diverge 
models. In Section 4, we solve stationary states and boundary fluxes for diverge models when 
vehicles have predefined routes. In Section 5, we discuss diverge models in various evacuation 
strategies. In Section 6, we demonstrate the validity of the proposed analytical framework with 
numerical examples. In Section 7, we summarize our findings and present some discussions. 

2 A multi-commodity continuous kinematic wave model of di- 
verging traffic flow and its Riemann problem 




-CXD x p = 



Figure 1 : An illustration of a diverge network 

We consider a diverge network with m > 2 downstream links and one upstream link, as shown 
in Figured] In this network, there are m+ 1 links and m paths. We differentiate all vehicles into 
P = m commodities according to their paths. We denote the link-commodity incidence variable by 
8p,a, which equals 1 if commodity p (p = 1 , • • • , m) uses link a (a= 1, • • • ,m + 1) and otherwise. 
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Thus P a = Lp=i &p,a is the number of commodities on link a: P a = 1 for a = 1 , • • • , m, and P a = m 
for a = 0. On a link a, the location is denoted by link coordinate x a G [X ai X a +L a ], where L a is the 
length of link a, and x a = X a and X a + L a are the upstream and downstream boundaries respectively. 
On the path of a commodity p, the location is denoted by commodity coordinate x p E [Xp, X p + LJ , 
where L p = £ fl Sp,aL a and we assume that there is no loop on a path, is the length of path p, and 
x p = X p and X p + L p are at the origin and destination respectively. If 8 p ^ a = 1 , we denote L pM as the 
distance from the origin of path p to the upstream boundary of link a, and x a and x p follows a one- 
to-one relation: if x p e [X p + L p ^X p + L p ]fl + L a ], then x p is on link a and x a = x p —Xp — L p ^ a +X a . 
That is, 8 p ^(x a — X a —x p +L pM +X p ) = for all a = 1, • • • ,m+ 1 and /? = 1, • • - ,m 

For commodity p, we denote density, speed, and flux by p p (x p ,t), v p (x p ,t), and q p (x p j) = 
pp(x p , t)v p (xp, t), respectively. From traffic conservation of commodity p, we can have the follow- 
ing continuous conservation equation 



dpp dqp 
dt dx„ 



0, 



(2) 



Haberman . 



1977 



Newelli 



1993). 



whose derivation is the same as that for single commodity (e.g. 
For link a, we denote density, speed, and flux by p a (x a ,t), v a (x a ,t), and q a (x a ,t) = p a (x ai t)v a (x ai t), 
respectively. Then we have that p a (x a ,t) = Lp8 p , a pp(xaj) and q a (x a ,t) = L p Sp, a qp(x ai t). Note 
that, p p (x a ,t) exists only when link a is on path p and p p (x a ,t) = p p (x p j) with x a = x p — 
L p ,a —X p + X a . It is the same for v p (x a j) and q p (x a ,t). We assume that traffic streams of 
different commodities on link a are homogeneous and share the same spe ed at the same loca 



tion and time. That is, we hav e the following speed-density relationships (|Greenshields . 



Del Castillo and BenitezL 



1995 



1935 



v p (x a ,t) = v a (x a ,t) =V(x a ,p a (x a ,t)). 



(3) 



Generally, V a (x a ,Pa) is non-increasing in p a , and Q(x a ,p a ) = p c y(x a: p a ) is unimodal in p a with 
its maximum as capacity at x a . We can see that conservation laws of multi-commodity flows in © 
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lead to the following LWR model 

^pa(Xa,t) + -^rPa(Xa,t)V(x a ,pa(Xa,t)) = 0, (4) 

which can work for inhomogeneous roads. Correspondingly, we can have the following traffic 
conservation equation for commodity p and x p G [X p + L p>a ,X p +L p>a +L a ] 

^Pp(xp,t) + -^-pp(xa,t)V(p a (x a ,t)) = 0, p=\,---,m (5) 

where x a = x p —X p — L pa +X a . For commodity p, the traffic stream evolves on the corresponding 
path, and we obtain a one-dimensional hyperbolic conservation law. However, all traffic streams 
interact with each other on the network, and we have a system of network hyperbolic conservation 
laws. We hereafter call © as a multi-commodity kinematic wave (MCKW) model of diverging 
traffic. 

We can see that traffic flow on a road network cannot be modeled by either one-dimensional 
or two-dimensional conservation laws, since vehicles of different commodities interact with each 
other on their shared links. In particular, for a diverge network with m downstream links, traffic 
streams of m commodities interact with each other on the upstream link. Traffic dynamics inside 
each link can be studied by the LWR models in (QQ) or ©, and the remained task is to study traffic 
dynamics at the diverging junction. Here we consider the Riemann problem for the MCKW model 
of diverging traffic in © with jump initial conditions. Without loss of generality, we assume 
that all links are homogeneous and infinitely long. For link a = l,---,m+l, we assume that its 
flow-density relation is q a = Q a (p a ), critical density p CjU , and its capacity C a . For the network in 
Figure [TJ we set X p — — °°, X p +L p = °o, and x p = is at the diverging junction for p = 1 , • • • , m; 
X a = and X a + L a = <*> for a = 1 , ■ • • , m; and Xq = — °° and X m+ \ + L m+ \ = 0. Therefore, L Pi o = 0, 
Lp,p — °°> an d Sp,a(x a —x p ) = for all a = and p = l,---,m. 

For commodity p = 1 , ■ • • , m, we have the following jump initial conditions: 

p p {x p ,0) = ^ . (6) 

P p ,r, x p e(0,+o<>) 
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Then upstream link and downstream link i = 1 , • • • , m have constant initial conditions: 



Po0o,o) = p = Y,Pp,L,x e (-°°,o), 



pi(xi,0) 



Pi = Pi,R, Xi E (0, +oo) , i=\, 



jn 



(7) 
(8) 



3 An analytical framework 



For link a — 0, • • • , m, we define the fo 



suppressed (lEngquist and Oshei . 



1980; 



lowing demand and supp 



Daganzo 



1995 



Lebacque . 



V func tions with all subscript a 



1996) 



D(p) = 2(min{p,p c }) 



Q(p), if p<p c 
C, if p > p c 

CP 



= [ P x(s)Q'(s)ds= f P max{Q'(s),0}ds 
Jo Jo 

<2(p), if p>p c 



(9) 



S(p) = 2(max{p,p c .}) = 

= C+ [ P (l-x(s))Q'(s)ds = C+ [ P min{Q'(s),0}ds, (10) 
Jo Jo 



C, if p < p c 

CP 



where %{p) equals 1 iff Q'(p) > and equals otherwise. 

Here we represent a traffic state in supply-demand space as U = (D, S) . This is different from 
many existing studies, in which traffic states are considered in p-q space. For the demand and 
supply functions in © and (TTOl) . we can see that D is non-decreasing with p and S non-increasing. 
Thus D < C, S < C, max{D,5} = C, and flow-rate q(U) = min{D,5}. In addition, D = S = C iff 
traffic is critical; D < S = C iff traffic is strictly under-critical (SUC); S < D = C iff traffic is strictly 
over-critical (SOC). Therefore, state U = (D,S) is under-critical (UC), iff S = C, or equivalently 
D < S; State U = (£>, S) is over-critical (OC), iff D = C, or equivalently S < D. 

In Figure Ob), we draw a supply-demand diagram for the two fundamental diagrams in Figure 
12 a). On the dashed branch of the supply-demand diagram, traffic is UC and U = {D,C) with 
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D < C; on the solid branch, traffic is OC and U = (C, 5) with S < C. Compared with the fun- 
damental diagram of a road section, the supply-demand diagram only considers its capacity C 
and criticality, but not other detailed characteristics such as critical density, jam density, or the 
shape of the fundamental diagram. That is, different fundamental diagrams can have the same 
demand-supply diagram, as long as they have the same capacity and are unimodal, and their criti- 
cal densities, jam densities, or shapes are not relevant. However, given a demand-supply diagram 
and its corresponding fundamental diagram, the points are one-to-one mapped. 




(a) 



(b) 



Figure 2: Fundamental diagrams and their corresponding supply-demand diagrams 

In supply-demand space, initial conditions in © and ([8]) are equivalent to (Here i = 1 , • • • , m if 
not otherwise mentioned) 



£/ (*o,0) = (A>, So), *oe(-oo,o), 
Ui{ Xi ,0) = (D h Si), Xi e(0,+oo). 



(11) 
(12) 
(13) 



In the solutions of the Riemann problem for © with initial conditions (|1 1H12I) . a shock wave or a 
rarefaction wave could initiate on a link from the diverging junction at x = 0, and traffic states on 
all links become stationary after a long time. We her eafter refer to these states as stationar y states. 



At the boundary, there can also exist interior states (Ivan Leeil . 



1984; 



Bult elle etal. . 



19981) . which 
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take infinitesimal space and only exist in one cell in numerical solutions. We denote the stationary 
states on upstream link and downstream link i by U Q ~ and £/ f + , respectively. We denote the 
interior states on links and i by Uo(0~,t) and Ui(0 + ,t), respectively. The structure of Riemann 
solutions on upstream and downstream links are shown in Figure |3j where arrows illustrate the 
directions of possible kinematic waves. Then the kinematic wave on upstream link is the solution 
of the corresponding LWR model with initial left and right conditions of Uq and U Q , respectively. 
Similarly, the kinematic wave on downstream link i is the solution of the corresponding LWR 
model with initial left and right conditions of and £/,-, respe ctively. 



Since vehicles' proportions travel forward along vehicles (|Lebac aue 



on the upstream link follow the First-In-First-Out (FIFO) principle (|PapageorgiouL 



1996), tra ffic dynamics 



1990). If the 



commodity proportions ^ are predefined and constant, as a result of the global FIFO principle, in 
the Riemann solutions we have 



(14) 



which serves as the First-In-First-Out principle (iPapageorgioul . ll99(J) . Also we have that, in the 
stationary state U Q , vehicles' proportions are the same as predefined ones. However, we could 
have different proportions in the interior state Uq(0~ ,t) and denote the corresponding proportion 
of commodity i by (0~ , t ) . 



Uo 



U - U o (0-,t) Ui(0+,t) u+ 

A L J L 



(a) 



(b) 



Figure 3: The structure of Riemann solutions: (a) Upstream link 0; (b) Downstream link i 



We denote #o-w as me A ux from link to link i for t > 0. The fluxes are determined by the 
stationary states: the out-flux of link is qo = q(U Q ~), and the in-flux of link i is qi = q(Uf). 
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Furthermore, from traffic conservation at a diverging junction, we have at stationary states 



q ^i = qi = q(Uf), q = q(U Q ) = ^q(U+). 

(=1 



3.1 Admissible stationary and interior states 



(15) 



As observed in (Hol den and Risebrol . 



1995 



Coclite et all 



20051) . the speed of a kinematic wave on 



an upstream link cannot be positive, and that on a downstream link cannot be negative. We have 
the following admissible conditions on stationary states. 

Theorem 3.1 (Admissible stationary states) For initial conditions in 4771) and ([72]), stationary 
states are admissible if and only if 



(D , C ) or (C ,S ), 



(16) 



where S Q < Dq , and 



(Q,Si) or (Df,Q) 



(17) 



where < Si . 



The proof is quite straightforward and omitted here. The regions of admissible upstream stationary 
states in both supply-demand and fundamental diagrams are shown in Figure 0] and the regions 
of admissible downstream stationary states are shown in Figure [51 From the figures, we can also 
determine the types and traveling directions of waves with given stationary and initial states on 
all links. In particular, the types of kinematic waves and the signs of the wave speeds can be 
determined in the supply-demand diagram, but the absolute values of the wave speeds have to be 
determined in the fundamental diagram. 

Remark 1. Uq = Uq and = Ui are always admissible. In this case, the stationary states are 
the same as the corresponding initial states, and there are no waves. 
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(c , A,) 



C„ D 




C a D 



(a) 



(b) 



Figure 4: Admissible stationary states for upstream link 0: marked by black dots 



m u, 



(a) 



£7,| X.~J?i'BL 



a, d o 



(b) 



Figure 5: Admissible stationary states for downstream link i: marked by black dots 
Remark2. Out-flux qo = min{D ,5 } < Dq and in-flux qi = min{D^ < St. That is. Do is 



the m aximum sending flow and 5/ is the maximum receiving flow in the sense of (iDaganzd . 



1995). 



1994, 



Remark 3. In (ILebacque and Khoshyaranl 



2005), a so-called "invariance principle" is proposed 
as follows: if Dq = Co, then q(U ~) < Dq; if Sf = Cu then q(U^) < Sq. We can see that Theorem 
13 . 1 1 is consistent with the "invariance principle". 



Corollary 3.2 For the upstream link 0, qo < Dq; qo < Dq if and only ifU = (Co, qo), and qo = Do 
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if and only ifU^ = (Dq, Cq). For the downstream link i, g, < Sa q v < Si if and only ifUf = (qi, Q), 
and qt = Si if and only ifU^ = (C,-, 5,-). That is, given out-fluxes and influxes, the stationary states 
can be uniquely determined. 

For interior states, the waves of the Riemann problem on link with left and right initial 
conditions of U Q ~ and Uq(0~ j) cannot have negative speeds. Similarly, the waves of the Riemann 
problem on link i with left and right initial conditions of C/ ; (0 + ,f) and cannot have positive 
speeds. Therefore, interior states Uq(0~ ,t) and C/;(0 + ,?) should satisfy the following admissible 
conditions. 

Theorem 3.3 (Admissible interior states) For asymptotic stationary states U ~ and Uf, interior 
states Uq(0~ ,t) and Ui(0 + ,t) in rf20l) are admissible if and only if 

{(Cn,5 n ) = L^rT, when U n is SOC; i.e., S„ < D^, = Cn 

o oo (lg) 

(Z)o(0~,?),So(0~,f)), when U Q is UC; i.e., D Q <Sq=Cq 

where So(0~ j) > Dq , and 

[ (Df.Q) = Uf, when t/+ is SUC; i.e., Df < Sf = Q 

Ui(0 + ,t) = { (19) 

[ (Di (0+ ,t),Si(0 + ,t)), when t/+ is OC; i. e., 5+ < D+ = Q 

where A(0 + ,t) >Sf . 

The proof is quite straightforward and omitted here. The regions of admissible upstream interior 
states in both supply-demand and fundamental diagrams are shown in ©, and the regions of ad- 
missible downstream interior states are shown in ©. From the figures, we can also determine 
the types and traveling directions of waves with given stationary and interior states on all links, 
but these waves are suppressed and cannot be observed, and we are only able to observe possible 
interior states in numerical solutions. 

Remark 1. Note that Uq{0~ ,f) = U Q ~ and L r ,(0 + , f ) = are always admissible. In this case, 
the interior states are the same as the stationary states. 
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Co D 



<>(/„- = u {o-,t) 



C a D 



(a) 



(b) 



Figure 6: Admissible interior states for upstream link 0: marked by black dots 




d D 



(a) 



(S+,C,) U,(<£,t) 




d D 



(b) 



Figure 7: Admissible interior states for downstream link i: marked by black dots 



Corollary 3.4 For upstream link 0, qo < Dq; q$ < Dq if and only ifUo(0~ , t) = U = (Co, qo), and 
qo = Dq if and only if Uq = (Do, Co), and Uq(0~ ,t) = (Dq(0~ ,t),So(0~ ,t)) with So(0~,t) > Dq. 
For the downstream link i, qi < St; qi < 5,- if and only ifUi(0 + ,t) = = (g ; -,C,), and qi = Si if 
and only ifU+ = (C,S ; ), and Uj(0 + ,t) = (A(0+,0,^(0 + ,0) with A(0 + ,0 > S/. 
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3.2 Entropy conditions consistent with discrete diverge models 

In order to uniquely determine the solutions of stationary states, we introduce a so-called entropy 
condition in interior states as follows: 

qt = F i (U (0-,t),U 1 (0 + ,t),---,U m (0 + ,t)^ 1 (0-,t),---^ m (0-,t)). (20) 

That is, the entropy condition uses "local" information in the sense that it determines boundary 
fluxes from interior states. In the discrete version of ©, the entropy condition is used to determine 
boundary fluxes from cells contingent to the diverging junction. Thus, 

Fi(u (o-, t),u 1 (o + ,t),---,u m (o + ,t)^i(o-,t),---Am(o-,t)) 

in (f20b can be consider ed as local, discrete flux functions. 



In (|Daganzo , 



1995), 



F(t/o(o-,0,^i(o + ,0.- } ^m(o + ,0^i(o-,0»--,&»(o-,0) 

was proposed to solve the following local optimization problem 



max {q } (21) 

UQ,u+,u (o-,t),Ui(o + ,t),-,u,n(o+,t)Mo-,t),-^n(o-,t) 



subject to 



qo<D (0 ,t), 
qi<Si(0 + ,t), 

£i(0~ ,t) = the proportion of vehicles choosing path i. 
Thus, we obtain the total flux as 

m S-(0 + t) 

F(tf (0V),^i(0V)r-^(0VU^ 

sal* ,t) 



In the literature, a number of other diverge models have been proposed. In (|Lebacqud . 11996), 



the upstream demand is split into commodity demands according to predefined turning proportions, 
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an d the in-flux of each downstream link is the minimum of its supply and commodity demand. 
In jjin and Zhangi LoCBb ). turning proportions were proposed to be d etermined by downstream 



1982T) . turning proportions were 



supplies when vehicles have no predefined routes. In (Sh effi et all 
proposed to be determined by downstream speeds in a myopic evacuation scheme. All these local, 
discrete diverge models can be considered as entropy conditions, so that we have corresponding 
continuous diverge models ©. 



3.3 Summary of the solution framework 

To solve the Riemann problem for © with the initial conditions in (fTTI)-(fT2l). we will first find 
stationary and interior states that satisfy the aforementioned entropy condition, admissible condi- 
tions, and traffic conservation equations. Then the kinematic wave on each link will be determined 
by the Riemann problem of the corresponding LWR model with initial and stationary states as 
initial conditions. Here we will only focus on solving the stationary states on all links, since the 
kinematic waves of the LWR model have been well studied in the literature. From all the condi- 
tions, we can see that the feasible domains of stationary and interior states are independent of the 
upstream supply, Si, and the downstream demand, D m+ \. That is, the same upstream demand and 
downstream supply will yield the same solutions of stationary and interior states. However, the 
upstream and downstream wave types and speeds on each can be related to 5, as shown in Figure 
HJd) and D m+ \ as shown in Figure |5£d). 



4 Diverge models with predefined turning proportions 

In this paper, we solve the Riemann problem for a diverging junctions with two downstream links; 
i.e., m = 2. In this section, we consider two entropy conditions, i.e., two diverge model. Here ve- 
hicles have predefined routes; i.e., are predefined constants, determined by vehicle route choice 
behaviors. We attempt to find the relationships between the boundary fluxes and the initial condi- 
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tions. 

qi = Pi(U ,U h U 2 ). (22) 

In contrast to local, discrete flux functions Fi(Uo(0~,t), C/i(0 - ,?), U2(0 + 1 t)), Pi(Uo,U\,Ui) can be 
considered as global, continuous. With the global, continuous fluxes, we can find stationary states 
from Corollary 13.21 With the solution framework in the preceding section, we can then find the 
kinematic waves of the Riemann problem of © with initial conditions (Uq, U\,U-i). 



4.1 Daganzo's diverge model 



In (|DaganzoLll995r) . a FIFO diverge model was proposed based on (|2T|) 



90 = min{D » (0 •' ) '|^7)'fe(F^) } ' <23) 

and a local FIFO principle 

qi = 6(0- Otfo- (24) 

Comparing (|24l) and (fl4l) . we obviously have <^(0~,?) = That is, the commodity proportions in 
the stationary state are the same as predefined. Thus in Riemann solutions, stationary and interior 
states have to satisfy (1231) , traffic conservation, and the corresponding admissible conditions. 

Theorem 4.1 For the Riemann problem of the MCKW model of merging traffic in (T5J) with ini- 
tial conditions in 4771) and rf72l) . boundary fluxes satisfying the entropy condition in ri2i1) . traffic 
conservation equations, and the corresponding admissible conditions are: 

q = mm{D ,^L,^}, (25) 
Si S2 

and qt = <^,go- The corresponding stationary and interior states are in the following: 

1. If Do > qo, the stationary state of the upstream link is SOC, and Uq(0~ , t) = U Q = (Co, go).' if 
Dq = qo, the stationary state of the upstream link is UC, Uq = (Do, Co), and Uq(0~ ,t) = Uq 
orUo(Q-,t) = (D (0-,t),S (0-,t)) with D (0~ \t) > D and S (0 - ,f) > A). 
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2. If 4> qo, the stationary state of downstream link i (i = 1, 2) is SUC, and Ui(0 + ,t) = = 
(qi,Q); if % — 1o> the stationary state of downstream link i is OC, = (Q,Si), and 
Ui(0+,t) = Uf orUi(0+,t) = (Di(0 + ,t),Si(0+,t)) with Z?i(0+ f) > S, andSi(0 + ,t) > S t . 

The proof of the theorem is given in Appendix A. The solutions of fluxes are illustrated in Figure 
[8l in which the starting point of an arrow represents the initial condition (Si,S2), and the ending 
point represents the solution (q\ , q 2 ) . That is, in region I, Dq < min,- %, and we have qi = ^Dq; in 
region II, < min{Z)o, %}, and we have q\ = S\ and q 2 = Si^-; m region III, < min{Z)o, 
and we have q 2 = S 2 and q\ = S 2 ^; on the boundary line between regions I and II, or the boundary 
line between regions I and III, = ^Dq; on the boundary line between regions II and III, q\ = Sj. 
We can see that, in region I, Dq < Si + S 2 , and q = min{D ,5i +52}. In regions II and III, 
qo < min{Do, Si + S2}. That is, due to vehicles' route choice behaviors, the capacity of the diverge, 
min{Do,Si +S2}, is generally under-utilized. 



Co . 

\ 




Si,qi 



Figure 8: The solutions of fluxes for a FIFO diverging junction 
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Comparing (1231) and (1251) . we can see that the global, continuous fluxes have the same func- 
tional form as the local, discrete fluxes. In this sense, the FIFO diverge model (1231) is "invariant". 
Hereafter, we consider a model invariant if and only if the global, continuous fluxes have the same 
functional form as the local, discrete fluxes. 



4.2 Lebacque's diverge model 



In (|Lebacqud . ll996l) . the following diverge model was proposed 

q t = imn{<^0V)A)((r,0A(0V)}, (26) 



and qo = q\ +q2- Compared with Daganzo's model (|23l) . this model is locally non-FIFO, and its 
solutions are the following. 

Theorem 4.2 For the Riemann problem of the MCKW model of merging traffic in (O with ini- 
tial conditions in 4771) and ( f72l) . boundary fluxes satisfying the entropy condition in ft26\l . traffic 
conservation equations, and the corresponding admissible conditions are the same as in rt25l) ; i.e., 

q = mm{D ,^,—}, 
SI S2 

and qi = ^iqo. The corresponding stationary and interior states are in the following : 

1. If Dq > qo, the stationary state of the upstream link is SOC, and Uq(0~ ,t) = Uq = (Co,qo); if 
Do = qo, the stationary state of the upstream link is UC, Uq = (Do, Co), and t/o(0~,?) = Uq 
orUo(0-,t) = (D (0-,t),S (0-,t)) withD (0-,t) > D and S (0~,f) > A)- 

2. 7/'| i > qo, the stationary state of downstream link i (i = 1,2) is SUC, and Ui(0 + j) = = 
(qi,Cj); if ^ = qo, the stationary state of downstream link i is OC, = (C,-,5 ; ), and 
Ui{0+,t) = Uf orUi(0+,t) = (A-(0V),S ; -(0V)) with Di(0+,t) > Si andS^O+j) > S t . 

3. The interior turning proportions £i (0~ , t ) can be determined by interior states and stationary 
states. 
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The proof of the theorem is in Appendix B. The solutions of boundary fluxes can also be illustrated 
by Figure [8] 

Compared with Daganzo's local FIFO model (|23T) . Lebacque's diverge model (|26|) does not 
satisfy the local optimization condition in (T2TI) or the local FIFO principle (l24l) and has interior 
commodity proportions different from the predefined ones. Due to the different functional forms 
of (|26l) and (|25l) . Lebacque's diverge model is not "invariant". However, Lebacque's model has 
exact the same fluxes, stationary states, and therefore kinematic waves as Daganzo's. That is, (l26l) 
converges to (l23l) asymptotically and continuously and yields globally optimal and FIFO solutions. 
In this sense, both models are "equivalent" globally and continuously. 



5 Diverge models for emergency evacuation 
5.1 A supply-proportional evacuation strategy 



We consider the a diverging rule proposed in (| Jin and ZhangL 12003b 



"1 r 



i, in which 



* = ^'■ 5,(0^)+^V.») }S ' (0+ -' ) - i=1 ' 2 ' <2?) 
In this diverging rule, vehicles do not have predefined routes and belong to the same commodity. 



This d iverge model was applied for emergency evacuation situations in a road network (lOiu and Jin , 



20081) . In this model, we have 

q = min{A)(0-,0,Si(0V)+S2(0V)}, 



and the turning proportions are time-dependent 

Si{0+,t) 



i = 1,2. (28) 



™ Si(0+ 0+52(0+ f) 

Thus in the Riemann solutions, stationary and interior states have to satisfy (l27l) . traffic conser- 
vation, and the corresponding admissible conditions. 
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Theorem 5.1 For the Riemann problem of the MCKW model of diverging traffic in (0) with initial 
conditions in 4771) and 4721) . stationary and interior states satisfying the entropy condition in < [27\) . 
traffic conservation equations, and the corresponding admissible conditions are the following: 



1. When Si + S 2 < Dq, Uf = Ui(0 + ,t) = (Q,Si) (i =1,2) and U = £/ (0~,f) = (C ,Si + S 2 ); 

2. When 5j +S 2 =D , Uf = U t (0 + ,t) = (Q,S f ) (i=l,2), U Q = (D ,C ), Uo(0~,t) = (D ,C ) 
or (D (0~,t),So(0~,t)) with Dq(Q~ ,t) > Dq andS (0~,t) > Dq when Dq < C ; 

3. When Si > <^A) (i = 1,2), U Q = U (0-,t) = (D ,C ), andUr = U t {0 + ,t) = (^D ,Q) 

4. When S\ + S 2 > Dq and Si < c ^_ c Dq (i : j = 1 or 2 and i ^ j), Uq = Uo(0~,t) = (Do, Co), 
Uf = (C u Si), Ui(0+,t) = (Cuj^Si), andUf = 17,(0+0 = (D -S„C ; ). 



The proof of the theorem is given in Appendix C. 

Corollary 5.2 For the Riemann problem of the MCKW model of diverging traffic in (0) with ini- 
tial conditions in 4771) and 4721) . boundary fluxes satisfying the entropy condition in H2 7D . traffic 
conservation equations, and the corresponding admissible conditions are the following: 



1. When S\+S 2 < Dq, qj = Si(i = 1 , 2) and qo = S\ + S 2 ; 



2. When Si > 



Q+c 2 



D (i=l,2), q { 



Ci+C 2 



Dq and qo = Dq; 



3. When S\ + S 2 > Dq and Si < 



Ci+C 2 



A) ( h j = 1 or 2 and i j), qj = Si, qj = Dq — Si, and 



qo = Dq. 



That is, for i, j = 1 or 2 and i ^ j, qQ = min{Do, S\ + S 2 }. 



min{5/, max{Do — Sj, 



Ci+C 2 



(29) 
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The solutions of fluxes in four different regions are shown in Figure HI in which the starting points 
of arrows represent the initial conditions in (Di,/^), and the ending points represent the solutions 
of fluxes (qi,q2)- In the figure, we can see four regimes: In regime I, both downstream links 
have OC stationary states; In regime II and IV, one downstream link has SUC and the other OC 
stationary states; In regime III, both downstream links have SUC stationary states. Comparing 
(1291) and (ITTT) . we can see that the evacuation diverge model is not "invariant". Compared with the 
diverge model in the preceding section, this model is optimal, since qo = min{Do,5i +S2}. 
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Figure 9: Solutions of fluxes for a supply-proportional emergency evacuation diverge model 



Corollary 5.3 IfUi(i= 1, 2) and Uq satisfy 



min{D/,5 ; } = min{S ; -,max{Do — S/, 
min{D ,5 } = min{5i +S 2 ,A)}, 



Dn 



C1+C2 



Q}}, 



then the unique stationary states are the same as the initial states, and traffic dynamics at the 
diverging junction are stationary. 
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5.2 A priority-based evacuation strategy 

Inspired by (|29| ), we propose a priority-based evacuation strategy (i, j = 1,2 and i ^ j) 

qi = min{5 ; (0+,0,max{D (0-,0-5 7 (0 + ,0,a ! -Z)o(0-,0}}, (30) 
where ce ; e [0, 1] and a\ + 0C2 = 1. 

Theorem 5.4 For the Riemann problem of the MCKW model of diverging traffic in (|5]) with ini- 
tial conditions in 4771) and 472]) . boundary fluxes satisfying the entropy condition in 001) . traffic 
conservation equations, and the corresponding admissible conditions are given by 

q t = min{S,-,max{D -S ; -,a/A)}}, ( 31 ) 

and qo = min{Do,5i +52}. 

The proof of the theorem is given in Appendix D. The solutions of fluxes (q\ , q?) from (Si , 52, Dq) 
are illustrated in Figure [TOl Clearly we can see that fluxes in (|29| ) can be considered as a special 
case when ^ = Q/(C\ +C2), and the priority-based evacuation diverge model is invariant. An 
extreme case is to give one downstream link an absolute priority for evacuation, e.g., <X\ = 1 and 
«2 = 0. This can happen when link 1 is shorter or less congestion prone. In this case the fluxes in 
(|3TT) become 

qi = min{5i,D }, 

q 2 = min{S 2 ,max{A)-Si}}. 

5.3 A partial evacuation strategy 

By a partial evacuation scenario, we mean that some vehicles have predefined routes and others do 
not. For example, £1 G [0, 1] and ^2 € [0, 1] are the predefined portions of vehicles choosing link 1 
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Figure 10: Solutions of fluxes for a priority-based evacuation diverge model 

and 2, respectively, but E,\ + ^2 ma y be smaller than 1. That is, the remaining portion 1 — £1 — ^2 
can take either route. For this scenario, we propose the following evacuation strategy (i, j = 1,2 
and i ^ j) 

qi = min{SK0V),i^j(0V) - Sj(0 + ,?),max{A)((r ,t) -Sj(0 + ,t), a,-D (0-,f) }}, (32) 

hi 

where a ; G 1 — £/] and oci + (X2 = 1 ■ 

Theorem 5.5 For Riemann problem of the MCKW model of diverging traffic in (|5]) with ini- 
tial conditions in 4771) and rf72l) . boundary fluxes satisfying the entropy condition in (TJ2|), traffic 
conservation equations, and the corresponding admissible conditions are given by 

q t = mm{Si,j-Sj-Sj,mait{Do-Sj,(XiDo}}, (33) 

and qo = q\ -\-q%. 

The proof of the theorem is omitted here. The solutions of (^1,^2) are illustrated in Figure [TT} in 
which there are six regimes. Furthermore, we can show that (1) qt > ^tqo', (2) When £1 + £2 = 1; 
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i.e., when all vehicles have predefined routes, (1331) is equivalent to (1251) ; (3) When <z,\ = £,2 = 0; 
i.e., when all vehicles have no predefined routes, ((331) is equivalent to (13"T1) . Therefore, (|32l) can be 
considered as a generalized diverge model, which encapsulates both normal and evacuation diverge 
model. In addition, the generalized diverge model is invariant. 




CiA) aiA, (l - 6)A> Ci D C 



Figure 1 1 : Solutions of fluxes for a generalized diverge model 



6 Numerical examples 



In this section, we numerically solve various diverge model and demonstrate the validity of our 
analytical results. Here, both links and 1 are two-lane mainline freeways with a corres ponding 



normalized maximum sensitivity fundamental diagram (IDel Castillo and Benitea . 

[0,2]) 



1995) is (p G 



G(P) 



p < 1 — exp 
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Link 2 is a one-lane off-ramp with a fundamental diagram as (p £ [0, 1]) 



,1,1 
1-exp -(--1, 



G(P) = ^{l-exp r _ V4Vp 

Note that here the free flow speed on the off -ramp is half of that on the mainline freeway, which is 
1. Thus we have the capacities Cq =C\ = AC2 = 0.3365 and the corresponding critical densities 
PcO = Pel = 2p c 2 = 0.4876. The length of all three links is the same as L = 10, and the simulation 
time duration is T = 360. Note that here all quantities in this section are normalized and therefore 
have no units. 

In the numerical examples, we discretize each link into M cells and divide the simulation time 
duration T into N steps. The time step At — T /N and the cell size Ax = L/M, with At = 0.9 Ax, 



satisfy the CFL condition (|Courant et all 



1928) 



At At 
Vf — = — = 0.9< 1. 
J Ax Ax 

Then we use the following finite difference equation for link £ = 0, 1,2: 



n n+l _ n n 1 At , n _ n \ 



where p" m is the average density in cell m of link i at time step n, and the boundary fluxes <fi m _\i2 
are determined by supply-demand methods. For example, for downstream links i =1 and 2, the 
out-fluxes are 



,M, 



where D" m is the demand of cell m on link i, S" m+l is the supply of cell m, and S" M+l is the supply 
of commodity i. For link 0, the in-fluxes are 

4,m-l/2 = ^{Dim-liSimh m=l,---,M, 

where Dq is the demand at the origin. Then the in-fluxes of the downstream links and the out-flux 
of the downstream link are determined by diverge models, which are discrete versions of (|20l ): 
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Qo,M+l/2 ~~ ?l,l/2 + #2,1/2 • 



We also track the commodity proportions in cell m of link 0, Q l m , as follows (|Jin and Zhang , 



2004 



:n+l _ Po,m %n , A? % m -l/2S>-l y0,m+l/2Sj, w 



Note that ^ is the predefined proportion of commodity i. 

In our numerical studies, we only consider Lebacque's diverge model (|26l ) and its invariant 
counterpart (|23l) . For Lebacque's diverge model, we have 

9j 1/2 = min{DS ;M ^,5^}, 

?0,M+l/2 = #1,1/2 +#2,1/2- 

In the invariant Daganzo's diverge model, we have (/, j = 1,2 and i ^ j) 

S n S n 

#0.M+l/2 = m M^O,M'J7r~'77T~}' 

n tin n 

#i,l/2 = ^iM%M+\ li- 
ft. \ Kinematic waves, stationary states, and interior states in Lebacque's 
diverge model 



In this subsection, we study numerical solutions of Lebacque's diverge model in (1261) . Initially, 
links and 1 carry OC flows with p\ = p3 = 1, and 30% of the vehicles on link diverge to link 2 
starting at t = 0; i.e., = 0.7, and ^2 = 0.3. The initial density on link 2 is p2 — 0. 1. That is, the ini- 
tial conditions in supply-demand space is U = Ui = (0.3365,0.2473) and U 7 _ (0.0500, 0841) . 



2004) 



Here we use the Neumann boundary condition in supply and demand (|Colella and Puckettl . 
Dq = Dq j , S'l M+l = S" M , and S 1 ^ M+l = ^2 m- Therefore, we have a Riemann problem here. 

In this case, < Do < Thus according to Theorem I4.2L we should have the follow- 
ing stationary and interior states U Q = Uo(0~ ,t) = (Cb,^|), = £/i (0 + , f) = {^S2,C\), and 
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= £/2(0 + ,f) = {C21S2) = (C2,C2). From the LWR model, there should be a back-traveling 
rarefaction wave on link connecting Uq to U Q , since So < %\ a forward-traveling shock wave on 
link 1 connecting to U\, since <S\\ and a forward- traveling rarefaction wave on link 2 
connecting Uj to U2. Furthermore, from (1351) . we should have that £1 (0~, f) = 0.5833. 




2 4 6 8 

x 2 



Figure 12: Solutions of Lebacque's diverge model (1261) : M = 160, iV = 6400. 



In Figure [121 the solutions of po, pi, and p2 are demonstrated with M = 160 and Af = 6400. 
From the figures, we can clearly see the predicted kinematic waves. In addition, we can ob- 
serve at t = T the approximate asymptotic values: Uq = Uo(0~ ,t) = (0.3365,0.2804), and p Q = 
po(0~v) = 0.8555; U+ = £/i(0+,f) = (0.1963,0.3365), and p+ = pi(0+,?) = 0.1963; and t/+ = 
U 2 {0 + ,t) = (0.0839,0.0841), and p+ = p 2 (0+,0 = 0.2436 w p 2c - These numbers are all very 
close to the theoretical values and get closer if we reduce Ax or increase T. That is, the results are 
consistent with theoretical results asymptotically. 

In Figure [T3J we demonstrate the evolution of the in-flux of link 1 and the proportion of com- 
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(a) Out-flux of link 1 (b) Proportion of commodity 1 in the last cell of link 
0.24 | , , , , 1 — | 0.7* , , , 1 1 — | 



cr 




Figure 13: Evolution of the out-flux and the density in the downstream cell of link 2 for Lebacque's 
diverge model (|26*1 ) 



modity 1 vehicles in the last cell of link for three different cell sizes. From Figure [TSTa) we 
can see that, initially, the out-flux of link 2 is minj^iDo, Si} = 0.2355, which is not the same but 
approaches the asymptotic in-flux ^2 = 0. 1963. Correspondingly the proportion of commodity 1 
vehicles in the last cell of link approaches the interior commodity proportion, as shown in Figure 
[T3f b). Similarly, as we decrease the cell size, the numerical results are closer to the theoretical ones 
at the same time. This figure shows that Lebacque's diverge model is not invariant, but approaches 
its invariant counterpart asymptotically. Note that the proportion of commodity 1 vehicles in any 
other cells of link remain constant at 0.7. 



6.2 Comparison of diverge models by Daganzo and Lebacque 

In this subsection, we compare the numerical solutions of Lebacque's diverge model (1261) with 
its invariant counterpart, Daganzo's diverge model (|23l) . Initially, links and 1 carry OC flows 
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with po = Pi = 1, and 30% of the vehicles on link diverge to link 2 starting at t = 0; i.e., 
(§i = 0.7, and ^2 = 0.3. The initial density on link 2 is p2 = 0.1. Different from the example in the 
preceding subsection, here we use the following boundary conditions: Dq = Dq j, S" M+l = S" M , 
and S'2 M+l = 0.05 + 0.03 s'm(nnAt/60). Thus we have a periodic supply on link 2. 




50 100 150 200 250 300 350 

t=nAt 

Figure 14: Difference in the solutions between Lebacque's diverge model (|26|) and its invariant 
counterpart (|23l 



We use p" m for the discrete density from Lebacque's diverge model (|26l) and p" m from its 
invariant counterpart (|23l) . Then we denote the difference between the two solutions by 

2 M 

e(nAt) = ZZKn-p'U^- (34) 

i=0 m= 1 

In Figure [H we can see that the difference decreases if we decreases the cell size. This clearly 
demonstrates that Lebacque's diverge model (l26l) converges to its invariant counterpart (|23|) . 
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7 Conclusion 



In this paper, we first introduced a continuous multi-commodity kinematic wave model for a di- 
verge network and defined its Riemann problem. Then, we introduced the supply-demand diagram 
of traffic flow and proposed a solution framework for the Riemann problem. In the Riemann solu- 
tions, each link has two new states: an interior state and a stationary state; and the kinematic waves 
on a link are determined by the initial state and the stationary state. We then derived admissible 
conditions for interior and stationary states and introduced entropy conditions consistent with var- 
ious discrete diverge models. In the analytical framework we proved that the stationary states and 
boundary fluxes exist and are unique for the Riemann problem for normal diverge, in which ve- 
hicles have predefined routes, and evacuation models, in which vehicles may not have predefined 
routes. With numerical examples, we demonstrated the validity of the solution framework devel- 
oped here and that Lebacque's diverge model converges to its invariant counterpart, Daganzo's 
diverge model, when we decrease the cell size. 

An important observation is that, for both (l26l) and (l27l) . fluxes computed by discrete supply- 
demand methods are different from the continuous fluxes. For example, the local fluxes from 
Lebacque's diverge model, (l26l) . are 

qi = min{^D ,5 ; }, i=l,2. 

When i^-Do > Sf, i.e., when the upstream demand is very heavy, we have qi = Sj. In this case, qi 
is not proportional to the turning proportion. Thus Lebacque's diverge model violates the FIFO 
principle. However, from the analysis in Section |4~21 and the numerical example in Section [6T2l we 
find that Lebacque's diverge model has the same continuous flux solutions as Daganzo's model, 
which observes the FIFO principle. Therefore, we conclude that Lebacque's diverge model is not 
strictly non-FIFO. As another example, for the supply -proportional evacuation model, at t = the 
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local fluxes from (l27l) are 



• r 1 ^0 l v 

min{1 '5 1 T^ }5i 



i =1,2, 



which are different from (1291) when only one downstream is SUC; i.e., when S\ + S2 > A) an d 
Si < q+C2^ ' H° wever ' the analytical results here suggest that the discrete fluxes converge to the 
continuous ones after a sufficient amount of time or at a given time but with decreasing period of 
a time interval. 

Comparing kinematic wave solutions of Daganzo's and Lebacque's diverge models, we find 
that, given the same initial conditions, they have the same stationary states and kinematic wave 
solutions, but different interior states. In this sense, interior states are essential to distinguish 
different diverge models. Numerical simulations in Section 16.21 also demonstrate the existence 
of interior states. Therefore, interior states are essential in understanding diverging traffic flow. 



This is different from the 



exist (|Jin and Zhang . 



2003a; 



, WR model for a homogeneous link, in which interior states could 



Jin et al 



2009) but are not essential to constructing kinematic wave 



solutions. 

Here we showed that both supply-proportional and priority-based diverge models can be con- 
sidered locally optimal evacuation strategies. But how to analyze kinematic waves arising in a 



speed-dependent evacuation model (|Sheffi et al 



19820 is subject to further investigations. In ad- 



dition to theoretical implications, this study, by improving our understanding of the formation and 
propagation of traffic congestion caused by diverging bottlenecks, could be helpful for developing, 
calibrating, and validating diverge models and associated emergency evacuation strategies in the 
future. For example, with different a,\ and 0,2, (1301 ) is a priority-based invariant diverge model, 
which can be used to evacuate vehicles to shorter or less congestion prone links without wasting 
the capacity of a diverging junction. In the future, we will also be interested in studying kinematic 
wave solutions of general junctions with multiple upstream and downstream junctions. 
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Appendix A: Proof of Theorem \4A 



Proof. From traffic conservation equations in (fT5l) . admissible conditions of stationary states, 
and the global FIFO principle (fT4)) . we have qi = ^qo < Si and qo < Dq. Thus, we have go < 
min{D ,|,|}. 
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Note that (1231) is equivalent to 

g = min{D (0 ,t), , }. 

SI S2 

We first prove ([25]). Otherwise, q < min{D , %}. (i) Since qo < D , from (OH) and COO), link 
is SOC, and Uo(0~ , t) = U Q = (C , qo) . (ii) Since q t = fyqo < Si, from dTT]) and £G1, link i is SUC, 
and Ui(0 + ,t) = Uf = (£#o,Q)- Hence from d23]) we have g = min{C , §-, < D < C . Thus 
g = min{|j-, %}. From the FIFO principle dHJ) we have = ^q = &(0~,£) min{^-, ^} > Q, 
and go = <?;/ > |S which contradicts qo < min{Do, %, 
We consider the following cases. 

(1) When only one of the three terms on the right hand-side of (1251) equals qo, we have (i) Do = 
qo < min{|-, (ii) |- = q < min{D , §}, or (iii) ^ = g < min{D , fj-}. Here we only 
show the solutions of stationary and interior states for (i), and solutions for (ii) and (iii) can be 
obtained in a similar fashion. When min ; j- > Dq = qo, from (|25T) we have qo = Do, and q t — 

Si 

$00 < From (OH) and (QJ), we have U Q = (D ,Cb), and £/ (0~,f) = (D (0~,0 ^(O - ,*)) 
with 5 (0~,0 > D . From Q7]) and (03, we have £7/(0 V) = £/+ = (&qo,Q)- Tnen from 
(|23l we have 

Q 

q = min{D (0",0,T i }- 

Since min,-^ > min,^ > Do = go, we have Do(0~,t) = go = Do, and f/o(0~,f) = Uq = 
(Do, Co). That is, the upstream and downstream interior states are the same as the corre- 
sponding stationary states. 

(2) When two of the three terms on the right hand-side of (|25|) equals qo, we have (i) Do = |j- = 
qo < (ii) Do = % = qo < fj-, or (iii) fj- = = qo < Do. Here we only show the solutions 
for (i), and solutions for (ii) and (iii) can be obtained in a similar fashion. When Do = < 

from Corollary!!! we have U Q = (D ,C ), and £/ (0 _ ,0 = (D ((T,?),So(0~,0) with 
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So(0- > D ; U+ = (Ci,Si), and Ui(0+,t) = (Z>i(0+ f),Si(0+0) with Di(0+ t) > Si; 
and t/ 2 (0 + ,?) = t/ 2 + = (g 2 ,C 2 ). Then from (f23]) we have 

• r . v Sl(0+,0 C 2 , 

go = min{A)(0 ,t), = ,-=-}, 

Si S2 

which leads to Do(0 - ,?) = Do or Si(0 + ,f) = ^\Dq = S\. In this case, we can have the follow- 
ing interior states for links 1 and 2: (a) Uo(0~,t) = U Q ~, and£/i(0+,?) = (Di(0+,?),Si(0 + ,?)) 
with Di (0+ t) > Si and Si (0+ ,t)>Sy, (b) U x (0+ , f ) = U+, and t/ (0~ , = (D ((T , , S ((T , 
with D (0", t) > D and S (0~,t) > D . 

(3) When all the three terms on the right hand-side of (T25T) equals qo, we have Do = |j- = & = go- 
From Corollary [3H we have t/" = (D ,C ), and £/ (0 _ ,r) = (D (0~,0,S (0~,0) with 
5o(0",0 > D ; f/+ = (C b Si), and U^O+j) = (Di(0+,?),Si(0+,0) with Di(0+ t) > Si; 
and £/+ = (C 2 ,S 2 ), and t/ 2 (0+ ?) = (D 2 (0+ ?),S 2 (0+,*)) withD 2 (0V) > S 2 . In this case, 
at least one of t/ (0~,?) = U Q , Ui(0 + ,t) = U^, and U 2 (0 + ,t) = £/ 2 + should be satisfied, and 
the other can be I7 (0~,f) = (D (0~,0,So(0 _ ,f)) with D (0 - ,f) > D and S (0~,t) > D , 
Ui(0 + ,t) = (Di(0+,z 1 ),Si(0 + ,f)) with Di(0+,f) > Si and Si(0+,*) > Si, or U 2 (0 + ,t) = 
(D 2 (0+,0,S 2 (0+,0) with D 2 (0+,f) > S 2 and S 2 (0+ ?) > S 2 . 



Appendix B: Proof of Theorem 33 



Proof. From traffic conservation equations in (1151) . admissible conditions of stationary states, 
and the global FIFO principle (fT4l) . we have g,- = <^go < S ; - and qo < Do- Thus, we have go < 
min{D ,|,|}. 

We first prove that qo is given by (|25l) . Otherwise, qo < min{Do, which leads to qo < Do 

and qi = fyqo < S ; -. From Corollary 13 .41 we have Uq(0~ ,t) = U Q = (Co, go) and t/,(0 + ,?) = = 
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(qi,Q). Then from (T26l) we have qi = min{^,(0~,?)Co,C,} < Sj < Sj. Thus qi = £,-(0~,f)Co, and 
<7o = Qi +12 = Co, which contradicts go < A) < Q). 
We consider the following cases. 



(1) When only one of the three terms on the right hand-side of (T25T) equals qo, we have (i) 
D = q < min{^, |}, (ii) | = q < mm{D , |}, or (iii) | = q < mm{D , |}. Here we 
only show the solutions of stationary and interior states for (i) and (ii), and solutions for (iii) 
can be obtained in a similar fashion. 

(i) When min, ^ > Dq = qo, from (T25l) we have qo = Do, and qt = fyqo < Sj. From Corol- 
lary m we have t/^ = (A),C ), £/o(0",f) = (Do(0-,0,So(0 _ ,0) with5 (0 _ ,0 > A), 
and f/,(0+,0 =Uf = (£#o,Q). Then from (US) we have 

9i = min{§{(0-,0/?o(0",0,Q} = ^(0",OOo(0~,0 <Q. 

Then we have go = A)(0~,f) = A)> and f/o(0~,f) = t/ ~ = (Do, Co). Further we have 
£;(0~,f) = ^j. In this case, the upstream and downstream interior states are the same as 
the corresponding stationary states. 

(ii) When g < mm{D , from <EU we have <? = § < A), 42 = S 2 , q\ = < s l- 
From Corollary [H we have U (0~,t) = Uq = (C , |), C/ 2 + = (C 2 ,S 2 ), t/ 2 (0+,*) = 
( J D 2 (0 + ,0,5 2 (0+,0) with D 2 (0+, f) >S 2 , and C/i(0 + ,?) = c/ i + = (f 5 2 ,Ci). Then from 
(|26l ) we have 

g 2 = min{^ 2 (0-,0C ,5 2 (0+,0}, 

9 1 = min{^(0-,OC ,Ci} = £i(0-,OCo<Ci. 

Thus we have 

€i(0-,0 = Irr, (35) 
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and 



qo = min{Co,S2(OV)+£i(OV)Co}=S2(OV)+£i(OV)Co<C . 
Thus, q 2 = S 2 (0+f) = S 2 , and U 2 (0+,t) = U+ = (C 2 ,S 2 ). 

(2) When two of the three terms on the right hand-side of (1251) equals qo, we have (i) Do = |j- = 
qo<-£, (ii) Do = ■£ = qo < or (iii) |^ = # = qo < Do. Here we only show the solutions 
for (i), and solutions for (ii) and (iii) can be obtained in a similar fashion. When Dq = -£< 
§, from Corollary[331 we have U Q = (D ,C ), and U (0~,t) = (D (0~ ,t),So(0~ ,t)) with 
5o(0",0 > £>o; ^ = (Ci,5i), and C/ 1 (0+,?) = (Di(0+,?),Si(0+,*)) with Di(0+ f) > Si; 
and t/ 2 (0+,?) = t/ 2 + = (q 2l C 2 ). Then from (O we have 

9l = min{^(0-,0A)(0-,0,5i(0 + ,0}, 

42 = min{fe(0",f)Do(0-,f),C 2 } = |2(0-,0^o(0",f) < C 2 . 

Thus Do(0~,f) = Do or Si(0 + ,f) = £iDo = Si. In this case, we can have the follow- 
ing interior states for links 1 and 2: (a) Uo{0~,t) = U Q , %j(0~,t) = §/, and Ui(0 + ,t) = 
(Di(0 + ,f),5i(0 + ,0) withDi(0+,f) > 5i andSi(0+,f) > Si; (b) £/i(0+,f) = Uf, &(0~,t) = 
&D /A)(0-,0» and £/ (0-,0 = (£>o(0- 0,S (0 _ ,0) withD (0-,?) > D and So(0~,0 > 
D . 

(3) When all the three terms on the right hand-side of (1231 ) equals qo, we have Dq = IX = || = 
From Corollary [3H we have £/ ~ = (D ,C ), and £/o((T,r) = (Do(0~,f),So(0~,f)) with 
5o(0",0 > D ; U+ = (Ci,Si), and C/i(0+,?) = (Z>i(0+ f),Si(0+0) with Di(0+ t) > Si; 
and U+ = (C 2 ,S 2 ), and U 2 (0+,t) = (D 2 (0+,t),S 2 (0+,t)) withD 2 (0+,*) > S 2 . In this case, 
at least one ofU (0~,t) = f/ ~, t/i(0 + ,?) = Uf, and t/ 2 (0 + ,0 = should be satisfied, and 
the other can be £/ ((T,0 = (D (0~,0,So(0 _ ,0) with D (0~,t) > D and So(0~,0 > A), 
£/i(0+,f) = (Di(0+,f),Si(0+,f)) with Di(0 + ,f) > Si and Si (0+0 > Si, or U 2 (0 + ,t) = 
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(D 2 (0 + ,0,S 2 (0V)) with D 2 (0 + ,t) > S 2 and S 2 (0 + ,t) > S 2 . Here &( 0_ :0 can be deter- 
mined once the interior states are determined. 



Appendix C: Proof of Theorem EL 



Proof. From traffic conservation equations in (1151) and admissible conditions of stationary states, 
we can see that 

qo < min{5i + S 2 ,D }. 

We first demonstrate that it is not possible that qo < min{5i + S 2 ,Dq} < min{Ci +C 2l Co}. Other- 
wise, from (flTT) and (fT9l) we have t/o(0~, t) = U Q = (Co, qo) with qo < Do; Since q(U^) +g(£/ 2 + ) = 
qo < Si + S 2 , then we have <?(£/+) < 5,- for at least one downstream link, e.g., q\ <S\. From (fT6l) 
and (fT8l) we have t/i(0 + ,?) = = {q\,C\). Then from the entropy condition in (f27l) we have 

go = min{Ci +S 2 (0 + ,t),C }, 

qi = min{1 'crrlko }Cl - 

Since qo < Co, from the first equation we have qo = C\ +S 2 (0 + ,t) < Co, and from the second 
equation we have q\ =C\, which contradicts q\ <S\. Therefore, 

qi +qi = qo = min{5i +S 2l D }. 

That is, the diverge model (1271) yields the optimal fluxes for any initial conditions. 

(1) When Si +S 2 < D , we have q = Si+S 2 < D . We have U (0~,t) = U Q = (Si+S 2 ,C ). 
Since q\ + q 2 = Si + S 2 and qi < Si, we have qi = Si, and = (C, Si). From (fT9b we have 
Ui(0+,t) = (Di(0+,t),Si(0 + ,t)) with A(0+,0 > S+ =S t . From ([27]) we have 

q = mm{Si(0 + ,t)+S 2 (0 + ,t),C } = S 1 +S 2 <Do<C , 

qi = ^ 1 ' ^H^ }5 ' (0+ ' f) = 5 ^ 
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Thus, 5/(0+, t) = Si < Si(0 + ,t). Then £/,-(0+,f) = U+ = (Q,Si). In this case, there are no 
interior states on all links. 

(2) When S\ + S 2 = Do, we have qo = Do, and qi = S ; -. We have U Q ~ = (Do, Co) and Uo(0~,t) = 

(D (0- , t) , D (0- , t) ) with So (0- , > D = D , and t/+ = (Q , Si) and Ui(0 + , t) = (S,-(0 + , t) ,5 ; (0+ , t)) 
with D, (0+, f) > Sf = Si. From we have 

q = mm{S 1 {0 + ,t)+S 2 (0 + ,t),D (0-,t)} = Si+S 2 = Do, 
We can have the following two scenarios. 

(2-i) If S 1 (0+,t) + S 2 (0+,t) > D (0-,0 = 5i +5 2 =D < S (0 _ ,0» then ^o(0",0 = U = 
(Do, Co) and there is no interior state on link 0. Moreover, we have 

51+52 S (0+ t) - s- 

which leads to Sj(0+,?) <Sj. From the assumption that5i(0 + ,?) +5 , 2(0 + ,f) >S{+S 2 , 
we have 5 ; (0 + ,f) = 5,-. Further we have f/,-(0+,?) = £7+ = (Q,Sj), and there are no 
interior states on links 1 or 2. 

(2-ii) If A)(0 _ ,f) > Si (0+0 +5 2 (0+,0 = Si +5 2 = D , S;(0+,0 = S t . Thus 17/(0+, = 
£/+ = (Ci,Si), and there are no interior states on links 1 or 2. Moreover, Uq(0~ ,t) 
satisfies So(0~,0 > Do and Do(0~,f) > Do- Thus there can be multiple interior states 
on link when Do < Co. 

(3,4) When Si+S 2 > D , then q = q l + q 2 = D . We have U Q = (D ,C ), and U (0~,t) = 
(Do(0~,0,Do(0~,0) with So(0~,0 > D = Do- For downstream links, at least one of the 
stationary states is SUC. Otherwise, from (flTT) we have E7+ = (C, Sf), and q\ +q 2 = S\ +S 2 > 
Dq, which is impossible. From (|27l) we have 



q = mm{S 1 (0 + ,t)+S 2 (0 + ,t),Do(0-,t)} = D <S l +S 2 , 
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If Si(0 + ,t)+S 2 (0 + ,t) < D (0-,t), then Si(0 + ,t) +S 2 (0 + ,t) =D < Si+S 2 and q t = 
Si(0 + , t) . This is not possible for the SUC stationary state = Ui(0 + ,t) = (qi, Q) with qi < 
Si < Q. Thus D (0~,0 < Si(0 + ,t)+S 2 (0 + ,t), D ((T, ?) = D < Si + S 2 , and £/ (0~,0 = 
U Q = (Do, Co). Hence for both downstream links 

Do 



Si(o+,o+s 2 (o+,o 



Si(o + ,0- 



(3) When S,- > j^Do (z = 1 , 2), stationary states on both links 1 and 2 are SUC with 
= £/ ; (0 + ,?) = (quCi) with < Sj. Otherwise, we assume that link 1 is SUC with 
E/i(0 + ,f) = C/+ = (^i,Ci) and link 2 is OC with t/+ = (S 2 ,C 2 ). Then 

* = c 1+ ^,,) 52( ° + '' ) -^ C2<52 - 



which is impossible. From (|27j> . we have 

M + C 2 

and/7 ; (0+,0 = t/+ = fe,C). 

(4) When Si + S 2 > S3 and S ; - < c ^ Do (z, j = 1 or 2 and i 7^ j), we can show that 
stationary states on links j and i are SUC and OC respectively with Uj~ = Uj(0 + ,t) = 
(qj,Cj) with qj < Sj, = (Q,S ( -), and Si(0 + ,t) > S/. Otherwise, Ui(0+,t) = Uf = 
(quQ) with 9; < S{, and 

D Q 

qi = c i+ Sj(^ t ) Ci -c^c 2 Do - Su 

which is impossible. Since at least one of the downstream links has SUC stationary 
state, the stationary states on links i and j are OC and SUC respectively. From (|27T) . we 
have a unique interior state on link i, t/,(0 + ,f) = (C/, D ^ s . Cj), and = Dq — S/. 



For the four cases, it is straightforward to show that (1291 ) always holds. 

44 



Appendix D: Proof of Theorem 15^4 



Proof. 

First (|3Q|) implies that 

q = min{Si(0V)+ 1 S 2 (0+, 0,A)(0",0}, 

which can be shown for three cases: (i) Si(0 + ,r) +52(0 + ,f) < Do(0~,f), (ii) S ; (0 + ,?) > ajSb(0 - ,?), 
and (iii) Si(Q + ,t)+S 2 (0 + ,t)>D (Q-,t) and 5,(0+, t) < onSo(0-,t). 

(1) When Si+S 2 < Dq, q = qy + q 2 < Si + S 2 < Dq < Cq. Thus the downstream stationary 
state is SOC with Uq = Uo(0~,t) = (Co, go)- In the following, we prove that q\ = Si, which 
is consistent with (TJTI). 



(i) Assuming that qi < Sj < Q, then the stationary state on link i is SUC with = 
Ui(0+,t) = {quQ). From d30j), we have 

qi = min{C, max{Co — Sj(0 + , t) , <x,Co}} = max{Co — Sj(0 + , t) , ce,-Co} < Q, 
qj = min{5j (0 + , t) , max{Co — Q , OCjCq } } ■ 

We show that the two equations have no solutions for either a ; Co < Sj(0 + : t) or OCjCo > 
Sj(0 + ,t). Thus qt = Si. 

(a) When ce/Co < S/(0 + ,t), we have a ; Co > Q — Sj(0 + ,r). From the first equation 
we have qi = a,Co. From the second equation we have qj = Sj(0 + 1 t) > OCjCo 
or qj = max{Co — C, CC/Co} > OCjCo. Thus qi + qj > Cq > Dq, which contradicts 
q < D . 

(b) When O/Cq > S/ (0 + , *) , we have a ; Co < Co - Sj (0 + , f) . From the first equation we 
have qj = Cq — Sj(0 + ,t). From the second equation we have qj = Sj(0 + ,t). Thus 
qi + qj = Cq, which contradicts qo < Dq. 
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(2) When Sj > CC/A), Dq — Sj < a ; A). In the following we show that qo = Dq and q t = a ; A), 
which is consistent with (f3Th . 

(i) If qo < Do, then the stationary state on link is SOC with U Q = Uo(0~ ,t) = (Co,<?o)- 
Also at least one of the downstream stationary states is SUC, since, otherwise, q\+q2 = 
S\ +S2 > Dq. Here we assume = C/,-(0 + ,f) = (q^Q). From (l30l) we have 

qi = min{C ; ,max{C -Sj(0 + ,t), a,C }} = max{C -Sj(0 + ,t), a,C }, 
qj = min{S/(0 + ,f),max{Co — Q, Of/Co}}. 

We show that the two equations have no solutions for either Co — Sj(0 + ,t) > (XjCo or 
Co-Sj(0 + ,t) < OiC . Thus q = D . 

(a) If Co — Sj(0 + ,f) > 05/Co, S/(0 + ,f) < 0/Go. From the first equation we have <ft = 
Cq — 5/(0 + ,f). From the second equation we have qj = Sj(0 + ,t). Thus qt + qj = 
Co, which contradicts qo < Do < Co. 

(b) If Co — S/(0 + ,f) < GCjCo, Sj(0 + ,t) > OyQ). From the first equation we have qj = 
OC,Co. From the second equation we have qj = Sj(0 + ,t) > (XjCo or qj = max{Co — 
Q, (XjCo) > (XjCo- Thus qi + qj > Co, which contradicts qo < Dq< Co. 

(ii) If q t < CtjDo < Sj < Q for any i = 1,2, then U+ = Ui(0+,t) = (quQ). From (EB we 
have 

qt = maK{D (0-,t)-Sj{0 + ,t),aiDo(0-,t)}<Q, 

qj = min{5j(0 + ,?),max{D (0 _ ,?) -C,a 7 -D (0 - ,?)}}. 

The first equation implies that (XiDo(0~ ,t) < oc,£>o; i.e., Do(0~ ,t) < Dq. In addition, 
D (0-,t) - Sj(0 + ,t) < (XiD . Thus, D (0 _ ,f) - C, < D - C, < D - ctjDo = (XjDo, 
and max{Z)o(0~,?) — C,, <XjDo(0~ ,t)} < Of/Do. From the second equation we have 
g ; < GCjDo- Thus g, + qj < Do, which contradicts qj + qj = Do- Thus g, > a,Do for 
i = 1,2. Since g,- + = Do, = a,-A). 
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(3) When 5,- + Sj > Dq and Si < CCjDo for i, j = 1 or 2 and i 7^ j. In the following we show that 
qo = Dq and = Sj, which is consistent with (T3TI) . 

(i) If qo < Dq, then the stationary state on link is SOC with U Q = Uq(0~ ,t) = (Co, go)- 
We first prove that at least one downstream stationary state is SUC and then that none 
of the downstream stationary states can be SUC. Therefore, q = Dq. 

(a) If none of the downstream stationary states are SUC, then q\ + q2 = S\ + S2 > Dq, 
which contradicts go < Dq. Thus, at least one of the downstream stationary states 
is SUC. 

(b) Assuming that qi < Sj, then Uf~ = Ui(0 + ,t) = (qj,Q). From (|3Q|) we have 

q t = min{Q, max{C - Sj(0 + , t) , aA}} = max{C - S/(0 + , t) , a,C } < S u 

which is not possible, since 5 ; < Of/A). Thus qi = Sj. 

(c) Assuming that qj < Sj, then U~j~ = Uj(0 + ,t) = (qj,Cj). Since qo — min{5/(0 + , t) + 
Cj,Co} < Dq, we have qQ = Si(0 + ,t) +Cj < Dq. From d30l) we have Si = qi < 
Si(0 + ,t). Thus Si + Sj < Si(0 + ,t)+Cj < D , which contradicts 5/ +Sj > Dq. 

(ii) If q t < Si, then Uf" = Ui(0 + ,t) = {quQ). From <f30]), we have 

qi = max{D (0~,t) -Sj(0 + ,t), aA)(0 - ,*)} < S; < a f A>, 
qj = min{5 7 (0 + ,0,max{D (0-,f)-C n a J D (0-,0}}. 

From the first equation we have that Do(0~,f) < Dq. We show that the two equa- 
tions have no solutions for either Do(0~,t) — Sj(0 + ,t) > 05,-Do(0~,f) or Dq(0~ ,t) — 
Sj(0 + ,t) < (XiDo(0-,t). Therefore q t = S t . 

(a) WhenD (0-,?)-Sy(0+,0 > aA)(0 _ ,0, we have 5/(0+ ?) < ajD (0~,t). Thus 
# = Do(0- - Sj(0 + ,t) and ^ = Sj(0+,t). Then 9j + qj = D ((T, ?) < Dq, 
which contradicts qi + qj = Dq. 
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(b) When D (O - , t ) - Sj (0+ , t ) < OfD (O - , t ), we have # = OfD (0 _ , J ) and D (0 _ , t ) - 
Q <Do(0~ ,t) -q t = (XjDo(0~ ,t). Thus^- < GCjD (0~,t), andqt+qj < D (0~ ,t) < 
Do, which contradicts q = D . 



48 



This figure "diverge20090226comparison.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/1009.4950vl 



This figure "lebacquediverge20090226contour.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/1009.4950vl 



This figure "lebacquediverge20090226interior.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/1009.4950vl 



